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ABSTRACT 

Globular clusters have proven to be essential to our understanding of many important 
astrophysical phenomena. Here we analyse spectroscopic observations of ten Halo 
globular clusters to determine their dark matter content, their tidal heating by the 
Galactic disc and halo, describe their metallicities and the likelihood that Newtonian 
dynamics explain their kinematics. We analyse a large number of members in all 
clusters, allowing us to address all these issues together, and we have included NGC 
288 and M30 to overlap with previous studies. We find that any flattening of the 
velocity dispersion profiles in the outer regions of our clusters can be explained by 
tidal heating. We also find that all our GCs have M/Ly < 5, therefore, we infer 
the observed dynamics do not require dark matter, or a modification of gravity. We 
suggest that the lack of tidal heating signatures in distant clusters indicates the Halo 
is not triaxial. The isothermal rotations of each cluster are measured, with M4 and 
NGC 288 exhibiting rotation at a level of 0.9 ± 0.1 km s" 1 and 0.25 ± 0.15 km s -1 , 
respectively. We also indirectly measure the tidal radius of NGC 6752, determining a 
more realistic figure for this cluster than current literature values. Lastly, an unresolved 
and intriguing puzzle is uncovered with regard to the cooling of the outer regions of 
all ten clusters. 
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1 INTRODUCTION 

Globular clusters (GCs) are often used as tracers of the 
gravitational potentials of galaxies and ga l axy clusters (e.g. 
Kissler-Patig et al J 1 19991: ICote et al.ll2003l; IWu fc Tremainei 
20061 : iQuercellini et all l2008i : iGebhardt fc Thomas! 120091 ). 
Although this has been applied to theoretical Galactic 
potentials (e.g. I Allen et all l200rj ). the actual Milky Way 
(MW) potential has not yet been analysed in this way. The 
tidal forces of spiral galaxies are thought to be strongest near 
the disc because the concentrated mass in that region (gas 
and stars) has a larger density gradient than the more slowly 
varying density of the dark matter (DM) halo. Interestingly, 
many distant MW objects such as GCs and dwarf galaxies 
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are known to be tidally stripped, despite being far enough 
from the Disc that they should not directly interact with 
it. For example, NGC 7492 is ~ 3kpc further from the 
Galactic centre than the most distant detection of the 
Monoceros Ring, an object on the very outskirts of the Disc 
(~ 22kpc: IConn et all l2007h . and exhibits clear evidence 
of tidal interaction with the Galaxy (|Lee et alj I2004T ) . In 
this paper we consider 10 GCs at varying Galactocentric 
and Planar distances, allowing the inference of properties of 
the potential of the MW for the first time, by looking for 
signatures of tidal heating in these clusters. 

Many globular clusters exhibit internal accelerations 
below ao ~ 1.2 x 10~ 10 ms ~ 2 , the level a t which either 
modified gravity (e.g. MOND : iMilgromll 19831 ) or dark matter 
is required to reconcile the observed kinematics of elliptical 
galaxies with theory. Near the tidal radius (rt) it is likely 
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that most stars in GCs feel accelerations below this level, 
making them an ideal tes t ing g round for low-acceleration 
gravity l)Sollima &: Nipoti l20ld . and references therein). 
Furthermore, if all GCs exhibit similar behaviour, Galactic 
influences can not be th e primary cause. In this final paper 
in the series [see also | Lane et al.l 120091 (hereafter Paper 
I) and lLane etafl 12010 J (hereafter Paper II)], we present 
the velocity dispersions and mass-to-light profiles of four 
GCs, namely M4, M12, NGC 288 (chosen for comparison 
with earlier studies) and NGC 6752, bringing the total 
for this project to 10. This sample allows statistically 
significant conclusions to be made on the dark matter 
content of Halo GCs, and on whether a modification of 
gravity is required to reconcile their internal kinematics with 
Newtonian gravitational theory. 

Our sample of GCs contains three close to the Galaxy 
(M55, M12 and M22), four at intermediate distances 
(NGC 6752, M4, M30 and 47 Tuc) and three that are distant 
(M68, NGC 288 and M53). We define 'close' to be R < 5kpc, 
'intermediate' as 5 < R < lOk pc, and 'distant' to be 
R > 10 kpc, following iHarrisI (|l996h . NGC 288 was chosen, in 
part, b ecause it is one of the GCs analysed bv lScarpa et al.l 
l|2007bl ) who found it to have a flat velocity dispersion 
profile, similar to that of Low Surface Brightness galaxies 
which are thought to be DM dominated through to their 
cores. Our targets were then analysed in separate studies 
(Papers I and II and the current paper) ensuring a mix of 
nearby, intermediate and distant GCs to ensure any Galactic 
influences, if any, would be clearly observed. See Table [1] for 
the estimated acceleration, due to the cluster, of the most 
distant cluster member for all ten clusters analysed in this 
project. Note that the three distant clusters all experience 
accelerations due to the Galaxy of ~ ao- 



2 DATA ACQUISITION AND REDUCTION 

We used AAOmega, a double-beam, multi-object spectro- 
graph on the 3.9m Anglo- Australian Telescope (AAT) at 
Siding Spring Observatory in New South Wales, Australia, 
to obtain the data for this survey. AAOmega covers a two- 
degree field of view, and is capable of obtaining spectra for 
392 individual objects over this field. We used 30 sky fibres 
used for optimal sky subtraction, and 5-8 fibres for guiding. 
The positional information for our t argets was taken from 
the 2MASS Point Source Catalogue jSkrutskie et alj|2006l ) 
which has an accuracy of ~ 0.1". 

Observations of M4 were performed on February 15-17, 
2008, with 1.5" - 2.5" seeing. The data for M12 were taken 
over two observing runs: 7 nights on August 12-18 2006, 
and a further 8 nights on August 30 - September 6 2007, 
both with mean seeing of ~ 1.5". NGC 288 was observed 
during the 2006 run and NGC 6752 during the 2007 run. For 
all observations we used the 2500V grating in the blue arm, 
resulting in spectra between 4800A and 5150A with A/AA = 
8000. In the red arm we used the 1700D grating, which is 
optimized for the Call IR triplet region. The red spectra 
cover 8350-8790A, with A/AA = 10000. This setup returns 
the highest spectral resolution available with AAOmega, and 
is suitable for measuring stellar radial velocities. We selected 
targets for this campaign by matching the J — K colour 
and K magnitude range of the red giant branch (RGB) of 



each cluster. To minimize scattered-light cross-talk between 
fibres, each configuration was limited to 3 magnitudes in 
range. 

We obtained 718, 2826, 1223 and 3664 spectra in the 
M4, M12, NGC 288 and NGC 6752 fields, respectively. 
Flat-field and arc-lamp exposures were used to ensure 
accurate data reduction and wavelength calibration. Data 
reduction was performed with the 2dfdr pipelin^E which 
was specifically developed for AAOmega data. We checked 
the efficacy of the pipeline with a comparison of individual 
stellar spectra. 

Radial velocities and atmospheric parameters were 
obtained through an iterative process, ta king the best 
fits to synthetic spectra from the library by iMunari et al.l 
(2005), degraded to the resolution of AAOmega, and cross- 
correlating this model with the observed spectra to calculate 
the radial velocity [a process very sim ilar to that used by the 
Radial Veloci t y Exp eriment (RAVE; Isteinmetz et aT]|2006l : 
IZwitter et al.l 12008) project]. We used the same spectral 
lib rary as the RAVE studies; this process is outlined in detail 
bv lKiss et all l|2007h . 



2.1 Cluster Membership 

We determined cluster membership using four parameters: 
the equivalent width of the calcium triplet lines, surface 
gravity, radial velocity and metallicity ([m/H]). Stars 
matching all criteria were judged to be members. Only 
stars having logg < 4.0 and log g < 4.6 were selected for 
NGC 6752 and M12, respectively, ensuring the majority of 
Galactic contaminants were removed before further selection 
criteria were applied. This probably removed some genuine 
cluster members but was necessary to ensure our sample was 
as free from Galactic field stars a s possible . 

For several clusters studied in lPaperU a cutoff of T e g > 
9000 K was necessary to remove hot horizontal branch (HB) 
stars. These have radial velocities with large uncertainties 
due to the calcium triplet in very hot stars being replaced 
by hydrogen Paschen lines. No cuts were made on T e s for any 
of the current clusters because no stars with T e « > 7425 K 
(for M4), T off > 5600 K (for M12), T cff > 7000 K (for NGC 
288) or T 6 ff > 5500 K (for NGC 6752) remained after our 
selection process. In total, 200, 242, 133 and 437 stars were 
found to be members of M4, M12, NGC 288 and NGC 6752, 
respectively. Figure [T] shows the relative locations of the 
observed stars and highlights those found to be members. 
Note that 19 stars in the M4 field were found to be members 
of the globular cluster NGC 6144. For M4, we find that 
88.0% of the selected members fall within within 2a of all 
selection parameters and 100% within 3o\ For M12 these 
values are: 94.2% and 100%, for NGC 288: 89.5% and 100% 
and for NGC 6752: 97.5% and 100%. Based on this we see 
no statistical reason to think there is significant Galactic 
contamination in our final samples. 



1 http:/ /www2. aao.gov.au/twiki/bin/view/Main/CookBook2dfdr 
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Figure 1. Distribution on the sky of the stars observed in the four fields, with axes in both degrees and parsecs from the cluster centre. 
Circled point s indicate stars that we determined to be cluster members (see text) . The large circle is the tidal radius of the cluster from 
lHarrisI II1996T ), with the smaller thinner circle in the lower right panel being our derived value for rt for NGC 6752 (see text). The large 
points in the upper right of the M4 field are the stars we determined to be members of NGC 6144. In each panel. North is up and East 
is to the left. 



3 RESULTS 

3.1 Tidal Radius of NGC 6752 

The tidal radius of NGC 6752 in the lHarrij \ 19961 ) catalogue 
is 55.'34 (shown as the large thick circle in the right panel of 
Figu repJI, however, this GC is kno wn to have a, collaps ed core 
(e.g. iRubenstein fc Bailvnl Il997l ) and lHarrisI jl996l ) warns 
against using tidal radii calculated from core parameters 
for such GCs. The most distant of our selected members 
throughout this project have, generally, been very close to 



the tidal radius. NGC 6752 is a clear outlier, with members 
only found t o within ~ 1/2 of the value of rt quoted by 
lHarrisI (|l996l ). Therefore, we propose an updated value of 
r t for NGC 6752 based on our membership selections. Note 
that M30 and 47 Tuc were not used to determine rt for 
NGC 6752 because of a paucity of stars observed in M30, 
and because stars were found out to the edge of the field of 
view of AAOmega for 47 Tuc. 

Based on the membership selections of the remaining 7 
clusters, the tidal radius is located 94.1±2.1% of the distance 



4 Richard R. Lane et al. 



to the most distant member from the cluster centre. Our 
most distant member for NGC 6752 is located at 33.87 pc, or 
29.' 11, therefore, rt 6752 = 27.'4±1.'7. The quoted uncertainty 
is based only on the standard deviation of the distance 
between the outermost member and the tidal radii of our 
clusters, so the true uncertainty in the value of r t is likely to 
be much greater than this. The large thin circle in the right 
hand panel of Figure [1] represents our derived value for rt- 
It should be stressed that this value is not a robust measure 
of rt but we sugg est it is a more realistic value than that in 



oi rt but w e suggest it is a n 
the lHarrisl (1 19961 ') catalogue. 



3.2 Metallicity 

1 



Our met allicity ([ Fe/H]) calibration method, discussed in 
detail in IPaper III , was used to determine the metallicities 
of each cluster. Briefly, the K magnitude of the Tip 
of the Red Giant Branch (-R'trgb) was subtracted from 
the K magnitudes of each star and plotted against the 
equivalent width of the calcium triplet lines to give a 
distance independent measure of luminosity. For M4, the 
^trgb value was taken from a J — K versus K Colour 
Magnitude Diagram based on 2MASS data within 5' of 
t he cl uster centre (i^TRGB = 5.3), for M12 from | Paus 
20061) (JTt rgb = 9.1), for NGC 288 from lDavidge~fe larri 



19971 ) and IValenti et all (120041) (K T rgb = 8.5), and for 

" ' 20041 ) (Ktbgb = 7.4). Linear 



NGC 6752 from Valenti etal 



fits to these data, combined with plotting [Fe /H] values vs 
T,W - AX for 47 Tuc and M55 ([Fe/H] from lHarrij|l996l . 
with A being the gradient of the slope above and X being 
K — Kthgb), allows a calibrator on [Fe/H] for other clusters. 

Figure [2] displays the robustness of this technique, 
with [Fe/H] values from this project plotted against those 
from the literature. Solid points are the clusters analysed 
in Papers I and II. From the current paper, the cross is 
M4, the square is M12, the triangle is NGC 288 and the 
diamond is NGC 6 752; the [Fe/H] literature values are 



.Johnson fc Pilachowskil (|2006l ). 
20001 ) and IZinnl (|l985 |) respectively . This 



from Kanatas et al.l l|l995l) 

IChen et al.l (|2000l ~ and fZmm qiaga i respectively , 
met hod is similar t o tha t outlined by I Cole et al] (|2004l ) 
and lWarren fe Cold (|2009l ). except we use the TRGB rather 
than the HB so it can be used for much more distant 
objects. A recent photometric study of 47 Tuc has revised 
the m etallicity of 47 Tuc to -0.83 (|Bergbusch fe Stetsonl 
2009). If this new value is adopted for our calibration, a 
maximum change in our calculated [Fe/H] values is -0.05 
(Kron 3), which is well within our uncertainty estimates. 
Calculated [Fe/H] values are also shown in Table [T] 



3.3 Rotation 

To measure the projected rotation of each cluster, we 
assumed an isothermal distribution. The rotations were 
measured by halving each by position angle (PA) and 
subtracting the mean stellar velocity of one half from the 
other. This was repeated in steps of 10° and the best-fitting 
sine function overplotted (Figure [3]). Note that for NGC 
6752 it was necessary to perform this in steps of 30° to 
avoid aliasing effects. The method results in an amplitude 
that is twice the projected rotation. Therefore, M4 exhibits 
rotation at 0.9 ± O.lkms" 1 , with an approximate axis of 
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Figure 2. [Fe/H] values derived from our method outlined in 
the text versus those from the literature. Solid points are those 
GCs from Papers I and II, and those from the current paper are: 
M4 (cross), M12 (square), NGC 288 (triangle) and NGC 6752 
(diam o nd). The [Fe/H] literature valu es are from iKanatas et all 
ll995h. Ijohnson fc Pilachowskil l|2006h . IChen et al] (|200Ct ) and 
IZinnl (119851 1 respectively. 



rotation of PA=70° - 250°, M12 at 0.15 ± 0.1 kms -1 , with 
an approximate axis of rotation of PA=40° — 220° (although 
this is effectively consistent with no rotation), NGC 288 at 
0.25 ± 0.15 kms -1 , with an approximate axis of rotation 
of PA=0° - 180°, and NGC 6752 shows no rotation to a 
level of 0.2 kms - . Our rotation measurem ent for M4 agrees 
well with that by IPeterson etaP (|l995h . who quoted an 
amplitude of 0.9 ± 0.4 kms -1 with an axis along the line 
PA=100° - 280°, although only for the inner 15' (~ r t /2, 
about 1/3 of the radius of our most distant member). 

For all clusters, we corrected the individual stellar 
velocity data for the measured rotation before calculating 
the velocity dispersions and M/Ly profiles. 



3.4 Velocity Dispersions 

The systemic velocities of each cluster were measu red using 
a Ma rkov Chain Monte Carlo (MCMC) method (jGregorvl 
2005), taking into account the individual velocity uncer- 
tainties on the stars, and providing the systemic velocities 
with associated uncertainties. A simple combination of the 
stellar velocities in each bin can provide a measure of 
systemic velocity and dispersion, although this does not 
take into account the individual velocity uncertainties. To 
fully incorporate these, we used a Bayesian MCMC-based 
analysis to provide a measure of realistic uncertainties of 
the velocity properties as a function of radius. Our systemic 
velocities agree very well with those from the literature (V r 
for all 10 clusters from this project are shown in Table 
[TJ, except for NGC 6144. Our survey did not sample this 
cluster well (19 stars were found to be members), however, 
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Figure 3. The rotation of each cluster calculated as the difference between the mean velocities on each side of the cluster along equal 
position angles, as described in the text. The best fitting sine function is overplotted, and a typical error bar is represented in the lower 
left of each panel. 



only 7 stars were analysed by iGeisler et alj l|l995h which 
may explain the discrepancy be tween the two values . The 
literature values are taken from ISuntzeff et alj (| 19931 ) (M4 
and M12).lRutledg . e et all (|l997h (NGC 288 and NGC 6752) 
and lGeisler et al.l (|l995h fNGC 6144). 

The velocity dispersions of our samples were calculated 
for annular bins centred on the cluster, each containing a 
similar number of stars (M4 « 20, M12 « 30, NGC 288 
« 20 and NGC 6752 w 40), centred on the cluster. The 
MCMC method described above was used to determine 
the dispersion in each bin, with the resulting velocity 
dispe rsion prof i les ov erplotted with the best-fitting line-of- 
sight IPlummerl |l91ll ) model: 



\R) = 



(1) 



Here, <to is the central velocity dispersion and r 3 is the 
Plummer scale radius (note that for Plummer models the 
unprojected half-mass radius is ~ 1.305 times the scale 
radius, however, for projected Plummer models such as in 
this paper r 3 is eq uivalent to the projected half-mass radius; 
iHaghi et all l2009h . The Plummer model is advantageous 
for our analysis because it is monotonically decreasing, so 
any flattening of the profiles would be discernible. It also 
allows for the calculation of the total mass of the cluster 
from the central velocity dispersion (<tq) and r s via (see 



iDeionghdl 19871 for a discussion of Plummer models and their 
application) : 



Mtot 



64agr s 
3ttG 



(2) 



Note that this model assumes the velocity distributions of 
the clusters are isotropic. How this assumption affects the 
overall conclusions of this paper regarding MOND and dark 
matter within GCs is far from obvious, and is a very complex 
pro blem which is beyond the scope of the current s tudy 
fsee lSpurzem et al]l20od : lGiersdl2006l : lKim et al.ll2008l . and 
references therein, for detailed discussions of the problem 
of isotropy in GCs). However, the anisotropies in each of 
the ten clusters analysed in Papers I, II and the current 
study are likely to vary greatly due to the large variation in 
Galactocentric distances and rotational velocities. Because 
all our clusters exhibit similar kinematic morphologies, with 
the possible exception of M4 (Section|4]), the anisotropy issue 
seems to have little impact on our results. 

We have chosen a Plummer model that does not include 
a tidal cutoff, over a more sophisticated Plummer model 
which includes a limiting tidal radius, because the tidal 
radii of many GCs are not well known (e.g. NGC 6752 
and 47 Tucanae; see Section ^. H and Paper II, respectively). 
Indeed, for many clusters we have found cluster members 
well outside literature values of rt- This means that our 
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Figure 4. Velocity dispersion profiles of each cluster from the current paper. The best fitting iPlummerl l ll91ll) model is overplotted and 
the derived scale radius, central dispersion and total mass is shown in each panel. 



model includes velocity dispersion information outside the 
tidal radii of our clusters. We have, therefore, not removed 
any velocity information for R> r t , which would reduce the 
accuracy of the model at large radii where a < ao- 

The velocity dispersion profiles, along with the 
total masses, scale radii and central velocity dispersions 
are presented in Figure U Except for M4, our mass 
estimates agree well w i th other studies (e.g. iMevlar] Il989l; 
iPrvor fc Mevlanl 1 19931 ; iKruiissen fc Mieskel 120091 ) , none of 
whom used Plummer models to calculate their estimates. 
For M4, we find a total mass about twice that of those 
studies [although this is reduced to a ~ 69% differenc e if the 
extremes of the uncertainties of IKruiissen fc Mieskel (|2009h 
and the current paper are taken]. Despite the Plummer 
profile fit being within the uncertainties, it is apparent that 
the outer five bins of M4 have nearly the same measured 
velocity dispersions. This increases the value of r s in the 
fitting of the profile, which is used to calculate the total 
mass (Equation [2J, which, in turn, leads to an inflated 
mass estimate. We attribute this apparent flattening of the 
dispersion profile to tidal heating (see Section HJ. 



Except for M4, none of the clusters discussed here have 
shown the appare nt flattening of the veloci ty dispersion 
profiles reported bv lScarpa et alj (2003, 2007a, b), indicating 
that neither a significant DM component, nor a modified 
theory of gravity, is required to explain their kinematic 
properties. This corroborates earlier results for 47 Tuc, 
M22, M30, M53, M55 and M68 i n Papers I and II 
and similar conclus ions are drawn by ISollima et al.l (l2009h 
for uj Centauri, b y Jordi et"al] l|2009f ) for Pal 14 and by 
ISollima fc Nipotil (|2010h for MOND theories in general. 
Through studies such as these, it is becoming increasingly 
apparent that neither DM, nor modified gravity theories, are 
necessary to explain the internal kinematics of GCs. 



3.5 Projected Mass-to-Light Profiles 

Dark matter causes larger stellar accelerations, and 
hence higher maximal stellar velocities, therefore, a good 
indication of whether a pressure-supported object like a GC 
is dark matter dominated is to measure its mass-to-light 
ratio. To calculate the M/Ly for our clusters we have used 
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Figure 5. Surface brightness data by iTrager et al . (1995) overplotted with best fitting Plummer profiles [note that these are fitted 
independently to the kinematic profiles in Section 13.41 and lead to a different value of r s (den oted r 3 j,), se e text]. The vertical lines 
represent the core radius (dashed), half-mass radius (dot-dashed) and tidal radius (dotted) from lHarris! 1 119961) . The solid vertical line in 
the lower right panel (NGC 6752) is the tidal radius we derived in Section [XT] 



the surface brightness data by (|Trager et al.ll 19951 . see Figure 
to which projected Plummer surface brightness profiles: 



I(R) 



Ltot 



r 2 
r s,L 



(r 2 



s.L 



+ R 2 ) 2 



(3) 



have been fitted (Figure [5j ■ These Plummer profiles were 
converted to solar lu minosities per sq uare parsec. Projected 
mass density profiles (|Peionghelll98^ 1: 



7T (r 2 + R 2 ) 2 



(4) 



in units of solar masses per square parsec, were then divided 
by the surface brightness profiles to produce radial mass- 
to-light profiles. The Plummer fits to the surface brightness 
data do not include tidal cutoffs, for the reasons discussed in 
Section \3. 41 Note that the kinematic and surface brightness 
models have been fitted independently, leading to two 



independent values of r s . The scale radius from the Plummer 
fits to the luminosity profiles (Equation [3] and Figure [5} is 
denoted r S: L', throughout this paper r s represents the value 
obtained from the cluster kinematics (Equation [JJ . These 
separate fits allow a radial, projected, M/Lv profile to be 
calculated: 



M/L v = ^ ( — 



r s ,L 



R 1 



+ R 2 



(5) 



The M/Lv profiles and mean values are shown in Figure 
[6] the thick line is the calculated M/Lv and the thin lines 
are their uncertainties. The M/Lv profiles at small radii 
deviate significantly from the mean. We interpret this as 
being due to the uncertainty in the measured luminosities 
and kinematics near the cores of the clusters. Crowding and 
confusion effects are inherent in luminosity and kinematic 
measurements of dense stellar fields, such as those near 
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Figure 6. Mass-to-light profiles the the four clusters analysed in the current paper. The thick line is the calculated M/Lv, the thin lines 
are their uncertainties and the vertical line is the value of r s . The quoted M/Lv value is only calculated for R > r s . None of these clusters 
have M/Lv 2> 1> furthering the argument that dark matter is not dominant. The deviation from a purely Plummer profile (Figure 
has the effect of decreasing the M/Lv of NGC 6752 by ~ 0.8 from its calculated value. 



the cores of GCs; it is clear in Figure [S] that there is a 
large spread of surface brightness measurements near the 
cores of all four clusters. Because of this uncertainty in 
core luminosities and kinematics, our mean M/Lv values 
were calculated for R > r s to ensure that these effects were 
removed. 

Dark-matter-dominated dynamical systems (e.g. ellip- 
tical and dwarf galaxies), exhibit high mass-to-light ratios 
(M/Lv > 10), whereas Ultra Compact Dwarfs, which 
have the same velocity d i spersion-luminosity rela tion as 
GCs (|Hasegan et all 120051 ; lEvstigneeva et al.1 120071 ), show 



no evidence for DM for M/Lv i$ 5- None of our clusters 
have M/Lv S> 1, therefore DM cannot dominate, and 
because none have M/Lv > 5 we see no need for any 
DM component. Similar conclusions were reached for all 
six clusters studied in Papers I and II. All results from 
this project indicate strongly that, in general, GCs do not 
contain large quantities of DM. We thoug ht it important to 
mention that iBaumgardt fc Mieskel (|2008l ) have shown that 
dynamically more evolved GCs exhibit lower M/Lv values, 
so the larger M/Lv of NGC 6752 should indicate that this 



cluster is dynamically "young". However, this is in direct 
contradiction with the c urrent understanding of NG C 6752 
having a collapsed core IjRubenstein fe Bailvnlll997l) . 

It is interesting to note that iFerraro et all (|2003l ) 
discussed the possibility of a large M/Lv value (~ 6 — 7) 
for the inner 0.08 pc of NGC 6752, because of the observed 
accelerations of millisecond pulsars near the core. Since we 
do not claim any knowledge of the M/Lv at those radii, 
this indeed remains a possibility, should the M/Lv increase 
further from r 3 toward the core. This intriguing possibility 
should be pursued by extracting more information on the 
kinematics of the core of NGC 6752, particularly in light 
of a newly recognised correlation that may b e useful for 
accur ately estimating the masses of GC cores l|Leigh et al.l 
l2009h . 



3.5.1 Correlation Between M/Lv and Luminosity 

47 Tuc is well known as having a bimodal distribution 
of various line strengths in both MS and giant stars (e.g. 
iNorris fc Freeman! [l979l ; lc annon et al.lll998l ; lHarbeck et al.l 
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Ta ble 1. All parameters derived from the kinematics of the clusters in this project i n order of decreasing metallicity (see text. [Paper j 
and lPaper III for literature comparisons), as well as the tidal radius from[Harris ( 1996|) and the estimated acceleration, due to the cluster, 
for the most distant cluster member. From left to right the columns are: cluster name, [Fe/H], systemic radial velocity, M/Lv, total 
cluster mass, rotational velocity, Plummer scale radius, central velocity dispersion, tidal radius and acceleration due to the cluster. Some 
parameters were not calculated due to low sampling, and 47 Tuc and M55 do not have calculated values for [Fe/H] because they were 
used as calibrators (see Paper II). V r , V I0 t and <ro are in kms -1 , Mtot is in 10 5 Mq, r s and rt are in pc and a is in ms~ 2 . The rt value 
for NGC 6752 in parentheses is that derived in Section |3. II 



Cluster 


[Fe/H] 


Vr 


M/L v 


Mtot 


Vvot 




CT0 


rt 


a 




47 Tuc 




-16.9 ±0.2 


4 6+°' 6 


nil 


2.2 ± 0.2 


7.8 ±0.9 


9.6 ±0.6 


56.1 


2.3 x 10" 


ii 


Kron 3 


-1.05 ± 0.10 




















M4 


-1.16 ±0.13 


71.5 ± 0.3 


1 6 +a5 
i - D -0.5 


9 9 + 0.7 
z - -0.6 


0.9 ± 0.1 


9.1 ±2.3 


3.9 ±0.7 


20.8 


3.8 x 10- 


11 


NGC 288 
Sgr 


-1.24 ±0.13 
-1.40 ±0.50 


-45.1 ±0.2 


2 1+ 1 ' 5 
^•■"■-1.3 


44+0-32 


0.25 ± 0.15 


3.7 ± 2.0 


2.7 ±0.8 


33.1 


6.2 x IO" 


12 


NGC 121 


-1.50 ±0.10 




















M12 
NGC 6752 


-1.50 ±0.13 
-1.62 ± 0.15 


-41.0 ±0.2 
-26.2 ±0.2 


i-i±S:S 
3-611:1 


5S+ ' 32 

2 +0 ' 6 
z ' u -0.6 


0.15 ±0.1 
nil 


1.5 ±0.8 
3.8 ± 1.1 


4.7 ±0.9 
5.7±0.7 


25.1 

64.4 (31.9 ± 2.0) 


1.1 x 10- 
2.7 x 10- 


11 
11 


M22 


-1.78 ±0.15 


-144.9 ±0.3 


4 7+ 1 ' 7 
'-1.7 


3.3+1? 


1.5 ± 0.4 


4.5 ± 1.5 


6.8 ±0.9 


27.0 


6.4 x 10- 


11 


M55 




174.8 ±0.4 


2 +0 ' 9 


1 4+0:5 


0.25 ± 0.09 


11.7 ±4.2 


2.7 ±0.5 


25.1 


1.9 x 10- 


11 


M53 


-1.99 ±0.10 


-62.8 ±0.3 


6.7+1'? 


5 2+ 15 
°- Z -1.4 


nil 


17.2 ±3.8 


4.4 ±0.9 


112.6 


5.8 x 10- 


12 


M68 


-2.06 ±0.15 


-94.9 ±0.3 


1 9+ i: ° 


"S7+0- 29 
u -°'-0.24 


0.6 ±0.4 


6.4 ± 2.0 


2.4 ±0.9 


90.0 


1.9 x 10- 


12 


M30 


-2.16 ±0.15 


-184.4 ±0.2 


i c+0.9 
•"••"-O.S 


90+ ' 51 
u - 9u -0.48 


nil 


2.3 ± 1.2 


5.0 ±0.9 


42.7 


8.1 x lO- 


12 


NGC 6144 




196.6 ±0.8 












82.2 







2003). The most popular explanation for this phenomenon 
is that the cluster has undergone multiple episodes of star 
formation. In Paper II we reported a rise in the velocity 
dispersion of 47 Tuc for R > rt/2, which we interpreted 
as evidence for evaporation. Here we present an alternative 
scenario for this signature: a past merger, which explains 
both the rise in velocity dispersion and the bimodality in 
line strengths, as well as its anomalously high luminosity 
compared with its M/Ly (see Figure [7]). 

The clusters from this project exhibit a clear trend 
between M/Lv and total luminosity (Figure [7] a trend 
which is generally attributed to low -mass star depletion, e.g. 
iKruiissen fc Portegies Zwartll2009l . and references therein), 
with 47 Tuc being an apparent outlier. Since there is no 
reason to think that 47 Tuc has a large DM component, 
a simple explanation is that 47 Tuc merged with another 
object of similar metallicity in its past. This merger would 
have increased the total luminosity of the cluster without 
altering its M/Lv significantly and caused bimodalities in 
both line strengths and velocities. If there has not yet been 
enough time for the two populations to mix thoroughly, 
this may be the cause of the observed rise in velocity 
dis persion. A detailed analysis of this scenario was presented 




dis persion. A detailed .• 
by (Lane et all |2010bl ) 



Absolute V Magnitude 

Figure 7. M/Ly from th is study vs absolute V band magnitude 
(taken from lHarrislll996h for all 10 clusters. Note the trend to 
higher M/Lv with increased luminosity. The open circle is 47 
Tuc. 



4 EVIDENCE FOR TIDAL HEATING 

When attempting to determine the reality of MOND using 
GCs, it is generally agreed that GCs at large distances 
from the Galaxy are most useful because the external 
acce leration imparted by th e Galactic tidal field is below 
00 (|Sollima Sz Nipotil l20ld . and references therein). An 
important question to ask, then, is at what Galactocentric 
distance does the external field become "negligible"? 
Furthermore, does the shape of the Halo have an effect on 



the GC dynamics, whereby distant GCs in certain regions 
of the Halo are affected above this threshold? One clear 
theoretical prediction is that the external field should heat 
the external parts of the GCs, thus increasing the velocity 
dispersion, up to, and inc luding, tidal destruction (e.g. Pal 5; 
lOdenkirchen et"al 120011 ). This is especially true during Disc 
crossings and at perigalacticon where tidal shocks strongly 
affect the dynamics of the cluster for short periods. We 
present here an examination of our data as an analysis of 
the Galactic tidal field. 
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Distance from the Plane (kpc) 

Figure 8. The ratio between velocity dispersions at th e tidal 
radiu s and the core vs distance from the Plane (taken from I Harris! 
Il996t) for the 10 clusters. Note the trend toward higher values of 
[crt/ero] toward the Galactic plane. The open circle is 47 Tuc and 
the open triangle is M4. 



Figure [8] shows how the ratio between the velocity 
dispersions at the tidal radius and core varies with distance 
from the Plane. The closer to the Plane, the greater the 
tidal effects from the Galaxy, and the larger the ratio. It is 
clear that for Rz > 3 kpc the tidal effects of the Galaxy 
are essentially equivalent at all radii. From this we can 
infer that the DM Halo exerts the dominant tidal force for 
Rz > 3 kpc. Furthermore, because the four clusters beyond 
R = 5 kpc (in increasing distance these are M30, M68 , NGC 
288 a nd M53) are in different locations in the Halo (|Harrisl 
1996) and on orbits with vastly different orientations to 
the Halo (I Allen et al.l l2006h . this indicates the poss ibility 
of a non-triaxial DM Halo (see lPenarrubia et alj feoog. for a 
discussion of the effects of halo triaxiality on the dynamics 
of GCs). Note that this is not strong evidence for the shape 
of the dark Halo, however, it is worth mentioning. 

M4 has a fairly flat velocity dispersion profile in the 
outer regions (Figure [4]). This cluster is also very close to 
the Plane (600 pc; lHarris! Il996h and has the largest value 
of [<Tt/(To] of any of our clusters (Figure [8]). Furthermore, 
its 3D space velocity, with respect to the Local Standard of 
Rest (jDinescu et al.l 1999), indicates that it is continuously 
interacting with the Disc because its orbit is nearly Planar. 
This, combined with its low M/Lv, strongly indicates that 
tidal heating is the cause of the flattening of the velocity 
dispersion in the outskirts of M4 rather than a substantial 
DM component. 

Figure [8] seems to indicate that GCs are tidally shocked 
within ~ 3 kpc of the Plane, then cool down, and stay 
cool, beyond that distance. Howe ver, the orbital periods 
calculated by iDinescu" et all (|l999l ) are about an order of 
magni tude shorter than the relaxation times in the iHarrij 
( 1996) catalogue for all of our clusters. Since the outskirts 
of GCs are much less dense than the cores, two-body 



interactions could not cause the outer regions to cool in 
less than the relaxation time of the cluster. In fact, no 
mechanism known to the authors can account for this rapid 
cooling in the outskirts of our GCs. Based on the M/Ly 
and velocity dispersion profiles of all 10 clusters, a large DM 
component is very unlikely. It is also very unlikely that we 
are seeing a MONDian effect in these outer regions because 
of the Newtonian nature of the velocity dispersion profiles. 
Addressing this interesting problem is beyond the scope of 
this paper, and more work is required to solve this intriguing 
puzzle. 



5 CONCLUSIONS 

In the current paper we have analysed four GCs (M4, 
M12, NGC 288 and NGC 6752) to determine their velocity 
dispersion and M/Lv profiles, bringing the total to 10 for 
this project. We have included GCs that have external 
accelerations extending from above ao down to ao and 
we find no deviation from our Plummer models at lower 
accelerations. Therefore, we see no indication that DM, or 
a modified version of gravitational theory, is required to 
reconcile GC dynamics with Newtonian gravity. This adds 
to the growing body of evidence that GCs are DM-poor, 
and that our understanding of weak-field gravitation is not 
incomplete. Within the stated uncertainties, the dynamics 
of all the s e clus ters are well described by purely analytic 
iPlummen (|l91lh models, which indicates that Newtonian 
gravity adequately describes their velocity dispersions, and 
we observe no breakdown of Newtonian gravity at ao ~ 
1.2 x 10~ 10 ms~ 2 , as has been claimed in previous studies. 

Despite this, we see the intriguing possibility of an 
unknown cooling process in the outskirts of GCs; the 
external regions of our GCs seem to cool much faster 
following tidal Disc shocks than the relaxation time of the 
clusters. Because it is highly unlikely that a MONDian 
process, or a significant DM component, is the cause of 
this cooling (based on our velocity dispersion and M/Ly 
profiles), further work is required to solve this puzzle. 
Furthermore, the lack of tidal heating events in the distant 
clusters provides some indication that the dark Halo is not 
triaxial. 

The Plummer model was also used to determine the 
total mass, scale radius, and M/Ly profile for each cluster. 
We find that none of our clusters have M/Ly ~S> 1, 
further evidence that DM does not dominate. We have 
produced M/Ly profiles, rather than quoting a single value 
based on the central velocity dispersion and central surface 
brightness. This method is used because it describes the 
M/Ly of the entire cluster, rather than only its core. 
This is particularly important for post-core-collapsed GCs, 
where crowding and confusion effects introduce significant 
uncertainty into luminosity and kinematic measurements at 
small radii. Within the uncertainties, our estimated cluster 
masses all match those in the literature except for M4, which 
we calculate to have a total mass about twice that of the 
literature values. The reason for this discrepancy is that the 
tidally heated cluster has an increased velocity dispersion in 
its outer regions, flattening the Plummer fit, increasing the 
value of r s , and therefore, increasing the mass estimate. 

Another important result from this study is the 
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measured rotations of our clusters. Of the four clusters 
studied here, M4 and NGC 288 show clear rotation, M12 
may have some rotation, and NGC 6752 displays no rotation 
signature. 

Throughout this project we have found similar results 
for the dark matter content, and Newtonian kinematics, 
of our 10 GCs, all at varying distances from the Galactic 
centre and Disc, including three that experience external 
accelerations due to the Galaxy of ~ oo- All data 
were acquired using the same instrument (AAOmega on 
the Anglo- Australian Telescope), reduced using the same 
pipeline (2dfdr), and analysed in the same way. This 
homogeneous approach is vital to a large project such 
as this, to ensure all systematics are accounted for in a 
similar fashion. Because of all these factors, our results from 
the three papers are strongly indicative that the current 
picture of globular clusters being dark-matter poor, and 
with dynamics explained by standard Newtonian theory, is 
correct. 
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